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Using the Density Matrix Renormalization Group (DMRG), we calculate critical exponents for 
the one-dimensional Hubbard model with open boundary conditions with and without additional 
boundary potentials at both ends. A direct comparison with open boundary condition Bethe Ansatz 
calculations provides a good check for the DMRG calculations on large system sizes. On the other 
hand, the DMRG calculations provide an independent check of the predictions of Conformal Field 
Theory, which are needed to obtain the critical exponents from the Bethe Ansatz. From Bethe 
Ansatz we predict the behaviour of the 1/L-corrected mean value of the Friedel oscillations (for the 
density and the magnetization) and the characteristic wave vectors, and show numerically that these 
conjectures are fulfilled with and without boundary potentials. The quality of the numerical results 
allows us to determine, for the first time, the behaviour of the coefficients of the Friedel oscillations 
as a function of the the Hubbard interaction. 

PACS numbers: 71.10.Fd, 71.10.Pm, 71.27.+a 
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I. INTRODUCTION 

Recent Bethe Ansatz studies of the one-dimensional 
Hubbard model with open boundaries subject to bound- 
ary chemical potentials or magnetic fieldsETQ have opened 
new possibilities to appiy-the predictions of Boundary 
Conformal Field TheorjiJlZI for the asymptotics of corre- 
lation functions to quantum_iuipurity problems. As in 
the case of periodic systems,Ht3 the corresponding ma- 
trix elements cannot be computed directly but must be 
extracted from the scaling behaviour of the low-lying ex- 
cited states. For generic filling and magnetization these 
finite size spectra allow the identification of the contribu- 
tions from two massless bosonic sectors associated with 
the spin and charge excitations in a Tomonaga-Luttinger 
liquid. 

A crucial step in these studies of systems with open 
boundary conditions is the correct interpretation of 
the finite-size spectra obtained from the Bethe Ansatz. 
These spectra determine both the bulk correlation func- 
tions, which are independent of the boundary fields, and 
the nonuniversal dependence of boundary phenomena 
such as the orthogonality exponent or X-ray edge sin- 
gularities on the strength of the scatterer. Hence both in 
the requirement of conformal invariance and in the anal- 
ysis of the finite-size spectra, the computation of corre- 
lation functions relies on assumptions which need to be 
verified by a more direct method, most notably by nu- 
merical calculations. Of course, comparison with exact 



results is also of interest for the numerical calculation: 
predictions of critical exponents for microscopic models 
allow the algorithms to be improved, which can then, in 
turn, be expected to produce better and more reliable 
results for more general systems. 

These considerations motivate our study of the Friedel 
oscillations for the single particle density and magnetiza- 
tion in the open Hubbard chain with boundary chemical 
potentials, which is described by the Hamiltonian 



L-l L 

n= - ( c L c -?+ i '°- + h - c -) + u ^2 n j,i n J,i 

",3=1 3=1 
— P ( n l,T + n l,l + n L,1 + n L,l) , (1) 

where the lattice has L sites, c'- a (cj- i<r ) creates (destroys) 

an electron on site j, rij t(7 — cj a Cj. a , and we have made 
the hopping integral dimensionless so that the Coulomb 
repulsion U and the on-site potential p are measured in 
dimensionless units. 

Numerical calculations of critical exponents in low di- 
mensional systems such as magnetic chains or electronic 
systems have a long history. Due to the need to consider 
systems of sufficient size, many earlier studies have used 
Quantum Monte Carlo methods to toat systems with pe- 
riodic boundary conditions (PBC)O More recently-the. 
Density Matrix Renormalization Group (DMRG)Bt3 
has become a new and powerful method especially suited 
to the study of one dimensional systems with open 
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boundary conditions (OBC).tjE-3 However systems, with 
PBC have also been studied with this method. Ejlla 

In the usual approach to the calculation of correla- 
tion functions with DMRG one considers large systems 
and averages over an ensemble of correlation functions lo- 
cated sufficiently far from the boundaries. In the thermo- 
dynamic limit, this procedure removes the Friedel oscil- 
lations due to the boundary and gives the bulk behaviour 
of the quantity in question. 

Here we want to make use of the existence of the ex- 
act solution of Eq. ([!]) in two ways in order to prepare 
the way for further extensions of the method: First, we 
use the quantities obtained from the Bethe Ansatz to 
provide checks of the numerical method at large system 
sizes. Second, we use the information contained in the 
oscillating behaviour to obtain more reliable results for 
the critical exponents. 

This paper is organized as follows. In Sec. [ll| we give a 
short description of the Bethe Ansatz and DMRG meth- 
ods, citing the relevant results from the Bethe Ansatz and 
Conformal Field Theory (CFT) concerning Friedel oscil- 
lations. In Sec. Ill, we study the Friedel oscillations of the 



density N(x) and the magnetization M(x). Combining 
the CFT results with those for noninteracting fermions, 
we obtain conjectures for the explicit form of the Friedel 
oscillations. After introducing the fit method used to 
obtain the exponents and coefficients from the DMRG 
results, we check the conjectures for two fixed densities 
and varying on-site interaction U. In addition, we study 
the dependence of the exponents and coefficients on the 
boundary potential p at one density. 



II. METHODS 

A. Bethe Ansatz 

The one-dimensional Hubbard model with OBC, Eq. 
(|l|), can be solved using the coordinate Bethe ansatz. 
The symmetrized Bethe Ansatz equations (BAE) deter- 
mining thepSpectrum of TL in the N e 
sector readEl □ 



N-t + iV r particle 



e^ 2L B c (k 3 )= JJ 



N„ 



Bs(Xa) J] 



A Q — sm(kj) + iu 



N 



sin(fcj) — \/3 + iu 
sin(fcj) — Xg — iu 



n 



A Q — Xg + 2iu 



-N, 



X a — sin(kj) — iu p N X a — X/g — 2iu 



B S (X) = 



X + 2iu 
X — 2iu 



(3) 



Since the BAE are already symmetrized and the solu- 
tions k — and A = have to be excluded, the energy of 
the corresponding eigenstate of Eq. (pi) is given by 



E 



1 — cos (kj). 



(4) 



In Refs. |2j-[4| the ground state and the low-lying excita- 
tions were studied for small boundary fields. In Ref. [Tj^ 
the existence of boundary states in the ground state for 
p > I was established. Bound states occur as additional 
complex solutions for the charge and spin rapidities. 

Here we will use the explicit form of the BAE, Eq. 
(||), to check the energy convergence of the DMRG re- 
sults for finite L. Furthermore, the expectation values of 
the density at the boundaries can be calculated from the 
derivative of the energy with respect to p (cf . Sec. [I B ) al- 
lowing another check of the numerics. Finally, the value 
of the magnetization at the boundaries for vanishing p 
can be calculated with a slightly modified Bethe-Ansatz 
(i.e. with a magnetic field at the boundary, see Refs. H,[|). 

Using standard procedures, the BAE for the ground 
state and low-lying excitations can be rewritten as lin- 
ear integral equations for the densities p c (k) and p s (X) 
of real quasi- momenta kj and spin rapidities A Q , respec- 
tively: 




K * 




(5) 



with the kernel K given by 



K = 





(i2u(X - sin(fc')) 



cos (k) a 2u (sin(fc) — A') 



-a 4u (X - A') 



(6) 



Here we have introduced a y (x) = y su +x i , and f * g 

denotes the convolution J^ A dyf(x — y)g{y) with bound- 
aries A = kg in the charge and A = Xq in the spin sector. 
The values of ko and Aq are fixed by the conditions 



1 



and 



-N R 



-Ni,...,^ (2) 



where we have defined u = £7/4 and identified the solu- 
tions k_j = —kj and A_ Q = — A Q in order to simplify the 
BAE. The boundary terms read 



ko 2 [N e - C c 
akpc = — 

-k L 



aA Ps — j , 

An L 



(7) 



(8) 



B c (k) 



1 — pe l 



sin(fc) 



sin(fe) 



where C c (C s ) denotes the number of-complex k (A)- 
solutions present in the ground state.liia In addition to 
the boundary terms in Eq. (||), the driving terms fP c and 
p s depend on whether or not the complex solutions are 
occupied or not. The explicit form can be found in Refs. 
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p|-[I|,|l9|. The presence of these 1/L corrections leads to 
the shifts 

9C p = \ I L I ° dk pc - 1 + 2C A ( 9 ) 



The correlation function G z aa ~ {M(x, t)M(Q, 0)) with 
magnetization M = iVf — iVj_ has the same critical be- 
haviour as G nn (x,t). Therefore (M(x)) has the same 
form as (N(x)), but with different coefficients A a . 




1 + 2C S 



(10) 



B. Density Matrix Renormalization Group 



where p c and p s denote the solution of Eq. (^|) without 
the 1/-7T driving term, i.e. the bulk system solution. 

Here we will be mainly interested in the exponents 
of the Friedel oscillations, given in Table § The quan- 
tity that determines, the critical exponents is the dressed 
charge matrix ZErEll 



Z C c Z c 



£cc(&o) &c(fco) 

£cs(Ao) £ss(Ao) 



(11) 



which is defined in terms of the integral equation 



i o 
o 1 



£,cc{k ) £,sc{k ) 
£ C5 (A') £„(A') 

(12) 



In Rcf. |2l] it was shown that the n-point correlation func- 
tions of the open boundary system are related to the 2n- 
point functions of the periodic boundary system. Thus 
the expectation value of the local density in the open 
system, (N(x)) can be extracted from the two-point 
correlation function (N(zi)N(z2)) v of the periodic sys- 
tem (see also Ref. where a spinless fermion model was 
considered). We can therefore use the results obtained 
in Refs. || and ^ for the density-density correlation func- 
tion. As a function of z c = x + iv c t and z s = x + iv s t 
(where v c and v s denote the Fermi velocities of the charge 
and spin sector, respectively), the asymptotic form of 
G nn (x,t) = {N{x,t)N{0,0)) is 



G nn (x, £) n e 
-Bj cos(2kp^x 



Bj cos(2fc Fa x + ^) 
\z c \^\z s \^i.* 
■Si) B n cos(2(fcir,t + kp,{)x + S n ) 



|27T, 



|27T,. 



|27n 



|27n, 



(13) 



with kF,i,i = 7rn ],i- The exponents are displayed in Ta- 
ble H Eq. (|l3|) shows the oscillating terms which are the 
most relevant ones asymptotically. For vanishing mag- 
netization the momenta fczrj. and /cf,t coincide, and one 
has to consider logarithmic corrections in x (see Ref. |23] ) 
- this case will notphe considered below. 

Following CardyEJ one has to replace \z c \ 2 — ► x and 
\z s \ 2 -> x to obtain (N(x)) from Eq. (jl|). The final 
result is 



(N(x)) =n e +A l 



■A, 



cos(2kp,iX + ipi) cos(2k f,i x + (p^) 
h 1 - A-> —. 

x 1l,c+~fi, s 1 x 7t,c+7t,s 

C0s(2(fc F | + fejP t |)x + (p n ) 



(14) 



y matrix renormalization group method 
3 ' has become one of the most powerful nu- 



The cLaa; 
(DMRGpl 

merical methods for calculating the low-energy proper- 
ties of one-dimensional strongly interacting quantum sys- 
tems. The expectation values of equal-time operators in 
the ground state, such as the local density or magnetiza- 
tion which interest us here, can be calculated with very 
good accuracy on quite large systems (on lattices of up to 
L = 400 sites in this paper). As we will see, access to such 
large system sizes is essential for the real space fitting 
method used to extract the coeffi cients and exponents of 
the Friedel oscillations (see Sec. |EII B[ ) . In the DMRG, 
open boundaries are also the most favourable type of 
boundary conditions numerically: for a given number of 
states kept (which corresponds to the amount of com- 
puter time needed) the accuracy in calculated quantitites 
such as the ground state energy is, in general, orders of 
magnitude better for open boundary conditions than for 
periodic boundary conditions oc3 

In this work the Unite system DMRG method is used: 
after the system is built up to a given size using a vari- 
ation of the infinite system method, a number of finite- 
system iterations are performed in which the overall size 
of the system (i.e. the superblock) is kept fixed, but 
part of the system (the system block) is built up. Opti- 
mal convergence is attained by increasing the number of 
states kept on each iteration, and the convergence of the 
exact diagonalization step is improved by keeping track 
of the basis transformations and using them to construct 
a good initial guess for the wavevectorEj For all calcula- 
tions shown in this paper, we have performed 5 iterations 
with a maximum of m = 600 states kept. The resulting 
discarded weight of the density matrix was 0(1O~ 7 ) and 
below. 

We illustrate the convergence of the algorithm explic- 
itly in Figs. 0-^. One finds that for all the parameters 
which are used in this paper the ground-state energy per 
site is accurate to 0(1O -6 ) or less, Fig. [j], while the ex- 
pectation value of the density at the first site (or at the 
last site, due to symmetry) is accurate to O(10 -5 ), as 
can be seen in Fig. |2[ 



TABLE I. Exponents 7 ajC and 7 QlS 
elements of the dressed charge matrix. 



as a function of the 
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FIG. 1. The difference between the ground state en- 
ergy per site calculated with the DMRG, Eg , and the exact 
BA-energy, Eq , as a function of the number of DMRG iter- 
ations for a typical system (L = 400, n e = 0.55, n-f = 0.35, 
U = 10 and p = 0). 




0123456789 10 
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FIG. 3. The difference between the expectation value of 
the magnetization at site 1 calculated with DMRG, Mf , and 
with BA, Mf, for n e = 0.55 and n e = 0.70, without boundary 
fields. 
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FIG. 2 
the density at site 



The difference between the expectation value of 
1 calculated with DMRG, Nf, and with 
BA, Nf, for the fillings n e = 0.55 and n e = 0.70 without 
boundary fields. 



The magnetization shows an analogous behaviour but 
with results correct up to O(10~ 4 ), Fig. ||. It is also 
interesting to note that for U < 10 there is no strong U- 
dependence of the accuracy of either the density or the 
magnetization expectation values. 



We now want to examine the effect of switching on the 
boundary fields simultaneously at the first and last sites 
on the quality of the DMRG results. In order to do this, 
we compare the mean density (Nx(p)) from DMRG calcu- 
lations with Bethe Ansatz results calculated in the ther- 
modynamic limit. Within the Bethe Ansatz, the mean 
density is calculated from the derivative of (Ti) with re- 
spect to the boundary field p. 



ST 1.0 




FIG. 4. The expectation value of {Ni) for different val- 
ues of the interaction U and boundary fields p on a L — 100 
lattice with n e = 0.55 and n-f = 0.35. The solid lines rep- 
resent the exact solutions in the thermodynamic limit, while 
the symbols are DMRG results. 



The numerical results, shown in Fig. [| on an L = 100 
lattice for electron density n e — 0.55, n-| = 0.35 and for 
two values of the interaction U , are again in very good 
agreement with the values for the thermodynamic limit. 
The difference Nf - Nf is now 0(1O~ 3 ). While this 
seems to be worse than the p = case, we have neglected 
finite-size corrections to (Ni) since we have compared to 
thermodynamic limit Bethe Ansatz calculations. If we 
explicitly take the 1/L corrections to the Bethe Ansatz 
values into account, we find agreement to 0(1O -5 ). This 
value is already smaller than the 1 /L 2 -corrections so that 
finally we state that (Ni) is accurate to 0(1O~ 4 ). 
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III. FRIEDEL OSCILLATIONS 



B. Fit procedure 



In general, the presence of an impurity or boundary in 
a one-dimensional fermion system leads to Friedel oscil- 
lations in the density, which have the general form 



8p(x) 



cos(2kpx + ip) 



(15) 



where the exponent 7 depends on the interaction. In ad- 
dition to numerical studies of these oscillations for spin- 
less fermional3 and Kondo-SystemsE£l, several theoreti- 
cal attempts have been made to clarify the role of in- 
teraction. Using bosonization it is possible to obtain the 
asymptotic exponents as a function of the interaction pa- 
ramctcrs. a i jidj corrections to the power-law behaviour of 
Eq. @.BE3 CFT results were used to calculate the in- 
teraction dependence of the exponent 7 for interacting 
spinlcss fermions. E3 Here we start with noninteracting 
fermions to obtain some conjectures for the connection 
between the explicit form of the Friedel oscillations and 
Bethe Ansatz results. These conjectures will then be 
checked using the DMRG results. 



A. Noninteracting Fermions 

By considering only spin-f electrons without any 
boundary potential, one can easily obtain the expecta- 
tion value of the electron density: 



N{x) 



iV T + i 



1 



sin ^Tra-tqpf- ) 



L + 1 2(L + 1) si / 



U+i 



(16) 



In the limit L — > 00 and x <C L the density is given by 



N(x) = n T 



» r 



n T - — 



sin ^27rx(ri| — " T L 2 )^ 
2nx 

27ra(n T - f ) 
2itx 



(17) 



with 6>g defined in Eq. (|). The density N(x) can also be 
calculated explicitly when the boundary field p = 1. It 
then has the same structure as Eq. (|l7|). 

If one assumes that the Friedel oscillations in the inter- 
acting system have an analogous structure to those in the 
non-interacting system, one can combine Eqs. ( |l4| ) and 
( |l7| ) to obtain the following conjectures for the finite-size 
shifts of the average density, average magnetization and 
the characteristic wave vectors in the interacting system: 



On 

-, m. 



(18) 



with 9 1 = 



6 A 

L 



L 



(19) 



= ®p~ ®pi ®n — 8% and 9 rn — d% — W s v 



Previously, several methods have been used to obtain 
asymptotic exnonents of correlation functions using nu- 
merical data.liilllj'La All of these methods use the in- 
dependence of the fourier-transformcd correlation func- 
tions near the relevant peaks k a (a —1,l,n) in Fourier 
space to extract the exponents. Due to the fact that 
only systems with periodic boundary conditions were 
considered, the k a were all independent of L. This L- 
independence seems to be crucial for these methods to 
work; we were not able to extract a reasonable exponent 
with any of these methods on a system with open bound- 
ary conditions. 

Therefore, we fit the DMRG results for N(x) and M(x) 
to the real-space test function 



/(*) = 



E 



A a sin(2fc Q a; + ip a ) 



ae{T,i,«} 



A a sin(2fc Q a; + ip a ) 
{L + 1 ~ x)^ 



(20) 



which explicitly includes the momenta as fit parameters. 
Here the second term is included due to symmetry. There 
are a total of 13 fit parameters in this function, a pro- 
hibitively large number to do a simultaneous fit of all pa- 
rameters. However, if we only consider systems in which 
the three peaks in the Fourier-spectrum are well seper- 
ated, there is an effective fit of 4 parameters to every 
peak. As we will see, the peak at k n is supressed for 
small U, reducing the number of fit parameters to 9 for 
only two peaks. The amplitudes A a will be assumed to 
be positive, with any sign given by the phase tp a . We fit 
to the magnetization M{x) and the density N(x) inde- 
pendently. 

The right side of Eq. (jll]) is only valid for x <C L. In 
addition, the CFT results are valid only asymptotically 
for large distances. As a consequence and compromise, 
we do not use the density information on the first five 
and last five lattice points. 

We perform the least squares fit in two stages. In the 
first stage, the start parameters of the subsequent fit are 
determined using simulated annealing. The final fit is 
performed using a combined Gauss-Newton and modi- 
fied Newton algorithm (using the NAG routine E04FCF). 
To estimate the fit error, 10 fits are performed for each 
system with 10% of the points randomly excluded from 
each fit. 
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C. Results 

Before discussing the results for the Friedel oscillations 
in detail, we make some general comments on the numer- 
ical results. As described in Sec. [II B, we calculate the 



quantities for the density and magnetization oscillations 
by applying a 13-parameter fit. Fitting to this many pa- 
rameters requires the use of large system sizes. While 
the numerical expense for the DMRG procedure grows 
linearly with the system size for a fixed number of states 
kept, the accuracy in the energy and in the local density 
and magnetization decreases with the system size, espe- 
cially in the Luttinger liquid regimeo We have compared 
the accuracy of the DMRG results with the accuracy and 
convergence of the fitting procedure for different lattice 
sizes from L = 300 to L — 500 and have decided that 
L = 400 yields optimal results for the amount of comput- 
ing power available. However, results within this range 
of system sizes are in agreement to within DMRG and 
fitting errors. 

Another important issue is the influence of the bound- 
ary potentials on the fitting method and on the Friedel 
oscillations (discussed in Sec. Ill C 3 ). As the boundary 
potential p isriacreased, bound states will develop at site 
1 and site L.ta In order to avoid these bound states in 
the fitting procedure, one has to enlarge the range in 
which the local density N(x) is disregarded from 5 (i.e. 
x, L — x — 1 < 5) at p = to about 20 at p = 9.9. 

The discussion of the next three sections will focus on 
comparing the BA/CFT predictions for the different fit 



parameters with the DMRG results, especially on check- 
ing the conjectures from Eqs. ( jl8| ) and (|l9|). We also 
compare the numerical results to the different exactly 
known values for different limits such as the limit of non- 
interacting fermions, t/ — ^ 0. 



1. Density n e = 0.55 

We first examine the Fourier transform of the local 
electron density N(x), defined as 



N(k, U) = J2 c os (k(x - h) N(x, U) 



(21) 



with k = -jj 1 - and j = 0, . . . , -j — 1. (Due to symmetry 



N(k,U) vanishes for odd multiples of -J.) The quantity 
N(k, U) is displayed in a three-dimensional plot in Fig. 
||(a). Distinct peaks at the three wave vectors, fcf, k^, 
and k n can clearly be seen. Note that we have chosen 
n e and nj so that these three peaks are well-separated. 
However, the peak at k n becomes very lightly weighted 
and therefore ill-defined for U < 1. In fact, we have 
found that it is not possible to locate the third momen- 
tum k n for U < 1 usi ng the 13-parameter fit procedure 
described in Sec. 1IIB. Therefore, we fit N(x) using only 
9 parameters for U — 0.1 and U = 0.5. We display 
the fourier-transformed magnetization M(k,U), defined 
analogously, in Fig. ||(b). Here the k n peak is even more 



N(k,U) 




(b) 



M(k,U) 




FIG. 5. Fourier transformation of (a) the density N(x) and (b) the magnetization M(x) for several values of U at the density 
n e = 0.55 and rif = 0.35. (The k = values are excluded.) 
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FIG. 6. Difference a = (a — a)L for the density (a = n) 
and the magnetization (a = m) as a function of U for 
n e = 0.55. The solid lines are the Bethe Ansatz conjectures 



poorly defined, and, in fact, is at best barely discernible, 
even at large U. Therefore, it is only possible to fit to two 
peaks using 9 parameters in the entire region from U = 
to U = 10. For these fit procedures, we have found that 
the mean-squared deviation a 2 = J2 X (-^( X ) — f( x )) 2 ^ s 
between O(10 -7 ) and 0(1(T 6 ) for the density and be- 
tween O(10 -6 ) and O(10 -5 ) for the magnetization for 
all U values. These limits will apply to all the fit results 
shown in this paper. 

In Fig. ^ we show the 1/L corrections of the mean val- 
ues n and m calculated with the DMRG and from Bethe 
Ansatz using the conjectures in Eq. (fl8|). One can see 
that there is quite good agreement between the two cal- 
culations for all U. 

The comparison of the exact asymptotic exponents at 
the different momenta with the numerical results is one 
of the most interesting and important features of this 
work because similar methods will then be able to be 
used to calculate properties not directly predictable with 
CFT/BA, and to treat systems that are not Bethe Ansatz 
solvable. 

The exponents at k± and fcj, extracted from the den- 
sity as well as from the magnetization data, are shown in 
Fig. 0. The difference between the fit-exponents and the 
CFT-prediction is less than 2% for the fits to the den- 
sity and less than 3% for the fits to the magnetization. 
As mentioned above, we can obtain an exponent for the 
peak at k n from the density fit for U > 1 only, due to its 
small weight especially for small U. At U = 1 the large 
error bars reflect a poor fit. Here we have only consid- 
ered U < 10, for which j n > 7|. A further increase in 
U would lead to a region where 7n < 7T- A crossover 
between these two regions will be seen for the n e = 0.70 
in the next section, for which it occurs at a somewhat 
lower U. 
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FIG. 7. Exponents for the three peaks as a function of U 
fitted for (a) the particle density and (b) the two peaks of the 
magnetization for n e = 0.55. The solid lines denote the CFT 
predictions from the Bethe Ansatz solution. 

Due to the fact that we obtain only three independent 
exponents from the fitting procedures, it is not possible to 
determine all of the elements of the dressed charge matrix 
Z. In fact, only the following combinations are relevant 
for the three exponents extracted: Z 2 C + Z 2 S , Z 2 S + Z 2 C 
and Z CC Z SC + Z cs Z ss r3 It would be possible to uniquely 
determine all of the independent elements of the dressed 
charge matrix with additional information from, for ex- 
ample, any susccptibilitycl or from another correlation 
function with a different set of critical exponents. Re- 
lationships between the elements of the dressed charge 
matrix Z and the parameters of the Tomonaga-Luttingcr 
model are given in Ref. ^9|. 

Within the framework of CFT, the amplitudes A a 
are completely undetermined. However, the form-factor 
approachEII may lead to explicit results for the amplitudes 
in the future. Fpr example, a conjecture of Lukyanov and 
ZamolodchikovEil concerning the amplitude of the spin- 
spin correlation functions of the XXZ_chain was recently 
confirmed by a fit to DMRG resultsO 

At this point, however, the fit results can only be com- 
pared to noninteracting fermions (U = 0), for which 
Ai = A | = j-. As can be seen in Fig. g, this value 
is in relatively good agreement (4% deviation) with the 
fit results. In addition, A n — for U = 0, in agreement 
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FIG. 8. Amplitudes A a fitted to (a) the density and (b) 
the magnetization as a function of U for n e = 0.55. The solid 
lines are guides to the eye. The arrow denotes the value for 
the noninteracting fermions (=^r)- Note that the amplitudes 
are, in general, different for the density and the magnetiza- 
tion. 

with the U = extrapolated value in Fig. pia). The large 
error bars in A n at U — 1 are due to the difficulty in fit- 
ting the k n peak for small U. Since the fitting procedure 
seems to work well for the exponents, and the amplitudes 
yield the correct [7 = limit, we feel that the calculation 
of the amplitudes is under good control. This is there- 
fore the first determination of the qualitative as well as 
quantitive behaviour of these amplitudes. 

The exact position of the momenta fcf and k± are fur- 
ther fit parameters. The l/L corrections to the ther- 
modynamic value are plotted in Fig. ^|. The fit values 
agree well with the Bethe Ansatz conjectures, except for 
the correction to k 



which deviates from the Bethe 



Ansatz value for U < 4. These deviations are probably 
due to problems with the fit. Note also that the Bethe 
Ansatz results for 6 a are correct only up to 0(1/ L). The 
momentum k n , which is not shown, is another indepen- 
dent fit parameter. The fit error in k n extracted from the 
fit to the density N(x) is rather large for U < 4. This 
is due to the fact that the peak at k n is not well-defined 
enough in this region to obtain the l/L corrections to this 
momentum. Nevertheless the agreement between the fit 
values and the Bethe Ansatz conjectures is very good for 




FIG. 9. Difference Q a ,8 = (n a V^)^ f° r a = li T as 

a function of U for the density n and magnetization m for 
n e = 0.55. The solid lines denote the Bethe Ansatz conjec- 
tures 6 a [see Eq. @]. 

U > 4. 



2. Density n e = 0.70 

In this section, we examine the same quantities as in 
the previous section at a density of n e = 0.70. A treat- 
ment of this density is interesting for a number of reasons. 
Since we use the same numerical parameters for both den- 
sities, we can examine the density dependence of the error 
in truncating the Hilbert space using the DMRG. This 
density is also interesting because the BA/CFT calcula- 
tions predict that the crossover between the 7„ and 7j 
exponents will take place within the range of interaction 
treated here, U = ... 10. 

The Fourier transforms of the local density N(x) and 
the magnetization Mix) are shown in Fig. Note that 
the momentum k n (wrapped back to the range to n) 
is now located between k^ and k-\. As can be seen in 
Fig. [H](a), the peak in N(x) at k n is not well-defined for 
U < 1, so the U = 0.1 and U = 0.5 data are fitted using 
9 parameters to fit two peaks. However, the k n peak in 
M(x) is now more well-defined than for n e — 0.55, as 
can be seen in Fig. |l^(b) , and it is now possible to fit to 
all three momenta for U > 2. For smaller values of U 
(U = 0.1, U = 0.5 and U = 1), a nine-parameter fit is 
again made to two peaks. 

The l/L corrections to the mean values of the density 
and the magnetization, shown in Fig. [ll], are in very good 
agreement with the BA conjectures, thereby providing a 
further confirmation of the predictions of Eqs. (18) and 
(|l9|). The exponents extracted from the fit are shown in 
Fig. O. The expected crossing of the two largest expo- 
nents at U f» 7 can clearly be seen. For 7^ and 71 ob- 
tained from the density fit, Fig. [121(a), the deviation from 



the CFT results is about 5% at most, with the largest er- 
rors occuring for U > 6, especially in 71 . As can be seen 
in Fig. |l^(a), the peak at fcj in N(k, U) gets weaker for 
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larger U, leading to a less effective fit. The agreement of 
the fitted exponents 7„ with the CFT predictions is much 
better, with a deviation from the BA/CFT values of less 
than 1% for U > 1. The exponents obtained from fits 
to the magnetization, Fig. [l^(b), show deviations of up 
to about 6% from the BA/CFT results. For U = 9 and 
U = 10, the deviation is largest and the exponents coinci- 
dentally take on the same value. Again, this is probably 
due to larger errors in the fit because the peak at fcj be- 
comes weaker at larger U. One can see that the peaks at 
fcf and k n are much less heavily weighted than the peak 
at fc| at large U. 

The amplitudes A[ and Af extracted from the fit to 
the density, shown in Fig. |l3|(a), decrease monotonically 
with increasing U . The U = values agree with the 
exactly known value of l/2w to within about 4%. The 
amplitude A n , on the other hand, increases with increas- 
ing U. Its U — > extrapolation agrees well with the 
value zero of the noninteracting fermions if the U = 1 
point, which cannot be very accurately determined, is ex- 
cluded. The behaviour and even the quantitative values 
of all three coefficients are quite similar to the n e = 0.55 
case shown previously in Fig. ||(a). The amplitudes ob- 
tained from the fit to the magnetization M(x) are shown 
in Fig. |l3](b). The amplitude A± behaves similarly to 
the n e = 0.55 case [Fig. ||(b)] in that it increases with 
increasing U, but Af shows different behaviour in that it 
reaches a maximum at U » 1 and then decreases. Both 
fit amplitudes yield the U — value of l/2ir to within 
6%. The amplitude A n for the summed momenta, which 



could not be determined for n e = 0.55, increases mono- 
tonically with U, and its U — > extrapolation agrees well 
with the value for noninteracting fermions, A n = 0. 

The \/L corrections to the momenta ki fitted for the 
density and the magnetization, shown in Fig. are in 
good agreement with the Bethe Ansatz conjecture. The 
agreement is also fairly good for 6^ n , although the error 
of the fit is rather large for small U. However, the 0\, m 
fit results do not match well with the conjecture. As we 
have seen in Fig. |l^(b), the k n and fcj peaks in M(x) 
have much lower amplitudes than the k-\ peak, leading to 
lower accuracy in the fitting procedure. The fit results 



a 




FIG. 11. Difference 6 a = (a — a)L for density (a = n) 
and magnetization (a = m) as a function of U for n e = 0.7. 
The solid lines are the Bethe Ansatz conjectures 8 a . 
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FIG. 12. Exponents for the three peaks of (a) the parti- 
cle density and (b) the magnetization as a function of U for 
n e = 0.7. The solid lines denote the CFT predictions from 
the Bethe Ansatz solution. 



for the exponents, Fig. |l2](b) also had a rather large de- 
viation from the CFT results in this regime. Thus it is 
not possible to confirm or deny the conjecture concerning 
the shift of fcj for M(x). 

For k n the situation is even worse. Both density and 
magnetization fits lead to large fit errors for U < 4. The 
deviation from the Bethc Ansatz conjectures is about 
0.08 in both fits, outside the range of the 1/L correction 
to 9 n . 

In summary, for n e = 0.70 the DMRG results for the 
mean density (magnetization, respectively) and the ex- 
ponents are in good agreement with exact results from 
BA/CFT, further confirming Conformal Field Theory. 

A detailed examination of the convergence of the 
DMRG shows that the numerical accuracy is actually 
slightly worse than the n e = 0.55 case, but this could 
be improved by increasing the number of states kept in 
the DMRG. We therefore expect to be able to apply 
these techniques reliably to obtain the boundary expo- 
nents and coeffients of other, non-Bethe- Ansatz solvable 
onc-dimcnsional models. 




0.00 — ' 1 1 1 1 1 1 1 1 1 

0123456789 10 

U 



FIG. 13. Amplitudes A a fitted to (a) the density and (b) 
the magnetization as a function of U for n e = 0.7. The solid 
lines are a guide to the eye. The arrow denotes the value ^ 
of the noninteracting fermions. Note that the amplitudes are, 
in general, different for the density and the magnetization. 
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FIG. 14. Difference 8 a ^ = (n a - for a =1, | as 

a function of U for the density n and magnetization m for 
n e — 0.7. The solid lines denote the Bethe Ansatz conjec- 
tures 6 a [see Eq. (B]. 



10 



ITP-UH-10/98 and WUE-ITP-98-015 



3. Effect of boundary potentials 

We now examine the effect of the boundary potential 
p on the Bethe Ansatz predictions. We set n e = 0.55, 
rij = 0.35, and U — 8 and consider six p values, for 
which four qualitatively different Bethe Ansatz solutions 
exist. For p = —0.5, 0, 0.5, no bound state is present in 
the Bethe Ansatz ground state; we examine both repul- 
sive and attractive fields. The values p = 2.6 and p = 6 
are in a region in which the Bethe Ansatz has two com- 
plex k solutions, one for each boundary potential at site 
x = 1 and x = L. The p = 6 ground state configuration 
contains, in addition, two complex A solutions. Finally, 
for p — 9.9 there are four complex k and two complex 
A solutions, corresponding to a bound pair of electrons 
at each end of the chain. Details of the structure of the 
ground state as a function of p are given in Ref. [l9| . 

The p-dependence of N(k,p) and M(k,p) at U = 8 
is not as strong as the [/-dependence of N(k,U) and 
M(k, U) found previously. Since we have chosen a fairly 
large U, the peaks in the density N(x) have enough 
weight to fit all three momenta. As was the case for 
p = 0, the peak at k n in the magnetization is not pro- 
nounced enough to be fitted at any p. 

The 1/L corrections to the mean values of the den- 
sity and magnetization, Fig. [l^, are again in very good 
agreement with the Bethe Ansatz conjectures, showing 
that Eqs. ( |l8| ) and ( [i"9| ) are valid even in the different 
physical regions described above. Within the BA/CFT 
calculations, the values of the exponents are indepen- 
dent of the boundary potentials p. This agrees with the 
DMRG results, which we do not show here: the range 
of the exponents varies by at most 2% from the exact 
values for p — (after a larger number of lattice points 
are discarded from the fits in order to avoid the bound 
states at the ends). 

The amplitudes also have no significant p-dependence. 
The density fit yields A l « A T m 0.045 and A n w 0.12, 
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FIG. 15. Difference 8 a = (a — a)L for density (a = n) and 
magnetization (a = m) as a function of p for U = 8, n e = 0.55 
and Tif = 0.35. The solid lines are the Bethe Ansatz conjec- 
tures for 6 a . 
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FIG. 16. Difference 6 a ,g = (n a — -^-)L for a =|,t,n 
as a function of p for the density n and magnetization m for 
U = 8, n e = 0.55 and n-j = 0.35. The solid lines denote the 
Bethe Ansatz conjectures 9 a [see Eq. (H)]. 

while the magnetization fit yields A± ~ 0.28 and Af w 
0.21. The absence of p-dependence at U = 8 suggests 
that the interaction dependence of the coefficients should 
be that of Fig. ||, independent of the boundary fields p. 

The effect of the boundary potential p on the shift of 
the positions of the peaks is much larger than the effect 
of varying U (compare Fig. [l6] with Fig. ^ . Since the 
fit results for all three k values agree very well with the 
Bethe Ansatz conjectures, the confirmation of Eq. ( |l9| ) is 
even more compelling than it was for the [/-dependence. 



IV. SUMMARY AND CONCLUSIONS 



We have carried out a detailed comparison between the 
exact Bethe Ansatz solution and Density Matrix Rcnor- 
malization Group calculations for the one-dimensional 
Hubbard model with open boundary conditions both 
with and without an additional chemical potential at 
both ends. A direct comparison of the ground state ener- 
gies as well as the density and magnetization at the ends 
of the chain has allowed us to estimate the accuracy of 
the DMRG on the large system sizes used in this work. 

We have then compared the behaviour of the Friedel 
oscillations in the local density and local magnetiza- 
tion calculated directly using the DMRG with Conformal 
Field Theory predictions for the asymptotic forms for 
which the exponents can be calculated using the Bethe 
Ansatz. We have performed this check for two different 
fillings, n e = 0.55 and n e = 0.70 for the case without 
boundary potentials, p — 0. We have obtained results 
consistent with the CFT predictions in all cases except 
those in which it is clear that the accuracy of the fitting 
procedure breaks down. Such a breakdown occurs when a 
particular peak in the fourier transform of the density or 
magnetization becomes lightly weighted and thus poorly 
defined. This occurs principally for the k n — (kp,-\ +&f,j.) 
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peak, especially at small U values. The good agreement 
between the CFT forms and BA values of the exponents 
and the DMRG calculations provides both a confirmation 
of the CFT predictions and a way to test the accuracy of 
the DMRG and of the effectiveness of fitting procedures 
for the Friedel oscillations. 

In addition, we have proposed a relation between the 
l/L corrected mean values in the density and magne- 
tization and the l/L corrections occuring in the BAE. 
This conjecture is supported by good agreement between 
mean values obtained from the fit to the DMRG data 
and the BAE results. 

We have been able to extract for the first time the in- 
teraction dependence of the amplitudes for the Friedel 
oscillations, a property not possible to calculate in the 
framework of the CFT, and have found the correct be- 
haviour in the U — ► limit. 

Finally, we have turned on boundary chemical poten- 
tials at n e = 0.55 and examined the p-dependence of the 
critical exponents, the amplitudes and our conjectures 
for the behaviour of the mean density and magnetiza- 
tion. The different p-regimes that we have considered 
yield qualitatively distinct Bethe Ansatz solutions that 
are physically connected with the formation of different 
types of bound states at the system boundaries. In agree- 



ment with BA/CFT predictions, we have found that the 
critical exponents are independent of p and that the influ- 
ence of p on the amplitudes is very weak. We have also 
found that our conjectures for the l/L corrected mean 
values of the density, magnetization and the wave- vector 
hold in all of the physically different p-regimes. 

The combination of analytical and numerical methods 
presented here has yielded new insights into both. The 
success of the numerical techniques will now allow the 
examination of more complicated systems that are not 
exactly solvable. Through comparison with the DMRG 
calculations, we have also been able to show that more 
information is contained in the BAE than is obtained 
from a direct interpretation via Conformal Field Theory. 
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